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The minimum rate of entropy production (MREP) and the least dissipation 
energy (LDE) principles are re-examined concerning continuous systems in sta- 
tionary nonequilibrium states. By means of simple considerations on coeffi- 
cients of phenomenological laws and by taking into account balance equations, 
an Onsager-like potential density for the stationary state is obtained. In the 
case of a single state variable it is found that the potential density describes a 
saddle-shaped surface around the stationary point in the domain of the force 
(state variable) and the flux. The saddle point may be considered a maximum 
along a LDE path (fixed force) and a minimum along a MREP path (fixed phe- 
nomenological coefficient). The latter allow for the Glansdorff-Prigoginc local 
potential. 



I. INTRODUCTION 

Thermodynamics establishes that, under certain conditions, systems in sta- 
tionary nonequilibrium states (SNES) have a minimum rate of entropy produc- 
tion (MREP) compatible with system constraints.^ The MREP principle (or 
criterion) requires that linear phenomenological relations between generalized 
fluxes and forces hold. This condition is satisfled when slight deviations from 
equilibrium take place uniformly in the whole system volume.^ Beyond this 
range some additional assumptions are required in order to apply a suitably 
modifled principle.'^ The matter is not so well defined in continuous systems 
when deviations from equilibrium are not uniformly distributed. Remaining 
near equilibrium, several problems can be assessed successfully provided the 
flux laws are suitably linearized.'' But in most case the linearization proce- 
dure presents some difficulties. Actually, according to statistical mechanics, 
linear laws hold when local thermodynamic equilibrium (LTE) conditions are 
satisfied.^ But, in this case linearity is meant in the sense of deviations from 
the local equilibrium. Thus, in general, phenomenological laws are not strictly 
linear. The matter is further complicated by a certain freedom in the choice 
of forces and conjugated fiuxes.^ Thus, for example, in the case of heat con- 
duction the MREP principle apparently works^ or fails^ according to the used 
conventions. 

Glansdorff and Prigogine''* formulated a general evolution criterion for ir- 
reversible processes which is independent of the actual form of the phenomeno- 
logical coefficients. Based on this criterion, they devised a procedure to define a 
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(local) potential for the SNES which coincides (unless a numeric factor) with the 
entropy production functional in the strictly linear range. Beyond this range, 
the local potential allows for a (generalized) MREP principle provided certain 
restrictions to the entropy production functional arc assumed.^ 

To give a thermodynamic framework to phenomenological laws, Onsager and 
Machlup proved a least dissipation energy (LDE) principle for near-equilibrium 
processes. ^"^'^^ In the assumption of Gaussian behaviour of extensive state vari- 
ables, the nonequilibrium processes are Markovian. A system displaced from 
equilibrium by a fluctuation decays back to equilibrium by following, on the 
average, the same empirical laws for the decay of the system displaced by an 
external constraint. According to this theory, the probability associated with 
a given decay path is to be searched for by means of an extremum principle 
of a suitable functional. The most probable path (with fixed forces) is the one 
established by phenomenological laws, i. e., by the LDE principle. 

The theory developed by Onsager and Machlup has given insights for the 
formulation of several Hamilton-based variational principles in the field of hydro- 
dynamics. We do not give accounts of the many works devoted to this topic (see 
for example ref. 9). In this paper we are interested in the possible connection 
between the Onsager (LDE) and the Glansdorff-Prigogine (extended MREP) 
principles for continuous systems. Specifically, we ask for a variational criterion 
including both the two principles to give a unified and concise representation 
of our knowledge about the SNES. To achieve its formulation we will pursue 
a simple task: the criterion must account for the (entropy and mass) balance 
equations regardless of the actual form of phenomenological coefiicients. 

To deal with this matter the paper is planned as follows. In Section 2 we will 
examine some features concerning the rate of entropy production as established 
by thermodynamics. In section 3 we will derive the variational criterion for 
the SNES. The extremum properties of a hmctional dependent on forces (state 
variables) and fluxes will be investigated. In the simplest case of a single (state) 
variable the density of the functional describes a saddle-shaped surface in the 
domain of force and flux. The extended MREP and LDE principles will be 
apparent by restrictions of the potential to special deviation paths. In the field 
of thermodynamics, restriction of the functional to the MREP path allows for 
the Glansdorff-Prigogine local potential. 



II Entropy production 

The definition of entropy production follows from the thermodynamic con- 
struction of the local entropy balance equation^'^ 

dps „ _ 

= + V • JsTOT (1) 

where s stands for the specific entropy , p for mass density , Jstot for total en- 
tropy fiux. The time derivative in the local frame is connected to the barycentre 
derivative by the relation 
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dt dt 

where v stands for the centre of mass velocity. By inserting eq. ^ in eq. H]) 
we define entropy flux 

Js = JsTOT - ps^r . (3) 

The starting point of the construction procedure is the Gibbs equation (assumed 
to hold in LTE conditions) in which energy and mass balance equations have 
been inserted. These latter read 

^ + v-a+p.v) = 5]i.,,e, (4) 

where pi stands for the density of the i-th system component, ji for diffusion 
flow, VijS,] is the rate of mass production of the i-th component in the j-th 
chemical reaction with reaction rate and stoichiometric coefficient Vij . The 
entropy flux can be written in general as 




where Jq stands for the heat flux and ipi for electrochemical potential of the i-th 
component^^. We are interested in those states where the barycentre derivatives 
vanish. Moreover, we consider the case of mechanical equilibrium where deriva- 
tives of V vanish everywhere in the system^. In these conditions, the density of 
entropy production is^^ 

i i 

where = ^^ij^^j stands for the affinity of the set of reactions involved 
in the production of the i-th component. When the centre of mass motion is 
determined by flows ji we have the condition Jj ~ 0. In this case the entropy 
production should be re- written in order to present only independent flows^^. 
Equation ([6]) can be used, for example, to investigate thermoelectric processes 
in semiconductors. In this case, chemical reactions are replaced by electron and 
hole excitation processes. 

Formally, as can be written as a product of generalized fluxes J and forces 

X 

k 

Fluxes are assumed to be functions of the forces (phenomenological laws): 



Jfc — E ^kl^l 



(8) 



By inversion of eqs. @ we have 



Xfc = MkiJi (9) 

where M ~ L^^. According to Curie's theorem, phenomenological coefficients 
are scalars in isotropic systems so that only fluxes and forces of the same vec- 
torial (or tensorial) character are related. The requirement of microscopic time 
reversibility of processes near equilibrium allow for symmetry of coefficients as- 
sociated with cross effects (this does not hold in cases of magnetic and Coriolis 
forces6>i2,i4) . 

= Lk^ (10) 

Owing to the detailed balance not all the possible definitions of forces and 
conjugated fluxes are allowed^ ( two important conventions consider "energy 
transport per molecule" or the "entropy transport per molecule" ^^). In this 
connection, we point out that a variational principle correctly formulated must 
not depend on the convention adopted. 

The entropy production of the whole system is (integrated over the system 
volume) 

fY^Jk-y^kd^'x (11) 

By considering the steady state as a reference state (Jfe ~ Jofc + (5Jfc,Xfc = 
Xofe -|-(5Xfc,), the first variation of P at the point {Xq, Jq) due to small deviations 
of forces is (constant Onsager's coefficients) 

S^'^P = / ^ Xofc ■ SJkd^'x + jY^ok- SXkd^'x (12) 

k k 

The second variation is 

fc 

Thus AP = S^^") P + S'^^'^ P ; no higher terms appear with the settings dealt with 
(see next section). Equation (fT3)) defines the excess of entropy production.^ If 
we consider equilibrium as the reference state {Xok ~ 0;Jofe = 0), we remain 

with AP = P = j Y.k^k- Xkd^x so that AP > and dAP/dt < (= at the 

stationary state) which explains the MREP principle. 

Let us we consider a simple example to clarify some points presented in the 
introduction. The diffusion under uniform temperature (for simplicity's sake 
we do not consider the Doufour and Soret effects) in absence of sources shows 
as = —j yp/p"' (we leave unspecified exponent n). A possible choice of the 
force and the conjugated flux is 



4 



X=--Vp 



(14) 



The variational calculus is applied to the functional P[p\ = J X{p)J{p)d^x by 
keeping fixed the variable at the system boundary: 



5P{p)^5j ±-{ypf<fx^ I 



— - V — 
dp dS/p J p^ 



ipSx = 



that is, 



d 
dg 



L 



2n- 



,2n+l 



ptl 



L 



Vp]}dp d^x = (15) 



Equation (|15p is apparently written in a complicated form, but it clearly shows 
that to account for the right transport equation (LVp/p" = const) we require 
that dL/dp = 0. The latter condition appears excessively restrictive since in 
some important cases (Fourier's law , Pick's law) Lap". A better choice may 
be: 

X -Vp J = j/p" = -L'Vp (16) 

where L' is a suitable constant coefficient. However, this choice does not match 
the transport features if the actual diffusion process requires V to be dependent 
on variable p. 



III. STATIONARY STATES BEYOND SMALL DEVIATIONS FROM 
EQUILIBRIUM. 

Henceforth, we will take SNES as the reference state from which deviations 
of variables are considered. To represent deviations of fluxes one can use expan- 
sions with respect to forces.^ But in this way information about phenomeno- 
logical coefficients are lost. In a more conservative representation deviations of 
coefficients are explicitly considered. Thus, by taking into account of expansions 

X, (p) = Xo, + 5(i)X, (p) + <5(2)x, (p) + ... (17) 

and 

(p) = Lozj + (p) + ^(^^Ly (p) + ... (18) 
we obtain from eq. ([5]): 

5(1) J. (p) = ^ Lo.,<5(i)X, (p) + ^ X« 5Uj ip) ■ (19) 

For convenience, in eqs (fT7|l and (fT8| we have indicated the functional depen- 
dence on state variables (for brevity symbolized by p) of all terms except the 
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quantities Xoj and Loij which pertain to the reference state. By taking into 
account that^ 

fj2^ok-S^'^Xk{p)d''x = (20) 

k 

eqs (|19|) and (fT2l) provide a way to represent SNES whatever conventions are 
used. Based on eq. ([^0]) we may write 

5(i)p[p] = J J2XokS^'^Jk {p)d'x (21) 

k 

which represents the variation of entropy production at the point po- In general, 
variation (S'-^-'P[p] is different from at the stationary state. To search for a 
potential showing vanishing variation at this state it is convenient to use the 
transformation 

/" ^ Xofe • <5(i) Jfe (p) d^x = ,5(1^ i /" ^ MonJk (p) ■ J; (p) d^x (22) 

•' k •' kl 

where Mq = L^^- Thus, by taking into account of eqs (HH), (|20p and (|22p we 
obtain 

-J^'^ - ^/ E^^o'«'Jfc(/^) • ^iiP)d'^'^ = (23) 

from which we define the functional 



Y.^k{p)-^k{p)d''x-- Y.Mokiik{p)-ii{p)d^x (24) 

A: •' kl 

and the associated density 

i>{p) ^Y.^k{p) • Xfe(p) - i ^ A%iJfe(p) ■ J((p) (25) 

According to definition ([2l]), we can rewrite eq. ([22]) as = (at SNES). 

It is worth to point out that the second deviation (5^^^ (p) does not present 
terms such as Xq^ • (^^^^ Jfc (p). This can be verified from eq. ((25)) . by taking 
into account eq. ^ written for the stationary state point, that is, 

S^^Hip) = 5]<5(^)Jfc(p) • <5(^)X,(p) - i5]Mofe,5(i)Jfe(p) •(5(i)j,(p) (26) 

k kl 

where the sums -^ofe ■ ^^^■'Jfe(p) multiplied by opposite signs are canceled 
out. In eq. (p6)) we have omitted the sum X]fc Jofc ■ '5'^'Xfc(p) since its integral 
vanishes as in eq. ([20]) . On these grounds, we can define the functional ^'[p, J] 
which has the form (|24p but in which the function 3k ip) is replaced by the 



6 



variable J^. Analogously, we follows the same rule to define the density ip{p,J) 
from eq. (|25p . Thus, the variational equation 

(5* = (27) 

be stationary at the SNES) can be solved with respect to the state variables 
(see eqs[T7] to[T9| or also with respect to the fluxes. In this case, the associated 
Eulero-Lagrange equations are 

{di^/d3i)sj,ES = ■ (28) 
which go back to the phenomenological laws at the SNES. The set of eqs (|28p 
devises a local extremum problem of the function -0 with respect to the fluxes 
by keeping the forces (state variables) fixed to their reference state values. To 
assure the existence of such an extremum it is sufficient for a general quadratic 
form /(r/) = J2ki ^okiVkVi to have a definite sign. In particular, ip {p,3) has a 
maximum (fixed forces) if the quadratic form is positive definite. 

In general, the functional does not show an extremum at SNES. However, 
as seen above, we can search for some special variation paths, crossing the 
SNES point in the domain of forces (state variables) and fluxes where this can 
be established in a restricted sense. To this end, we extend the problem (^5)) 
to any set of flxed forces lying within a small interval around the SNES point. 
The locus of maxima thus found satisfies the relation 

Xfe = ^A/ofciJi (29) 
I 

which by inversion gives 

Jfc=^iofciXi. (30) 

I 

Note that the points represented by eqs (|30p does not satisfy the LTE conditions 
which are fitted by eq. ([5|). Thus, deviations of fiuxes defined by 

^J^P) -I]iofei^Xi(p) (31) 
I 

does not introduce any additional condition such as SLf^i (p) — 0. In the view of 
a local LDE principle, deviations (pij) may be considered as fluctuations from 
the (SNES) local equilibrium. From substitution of eq. (|501) into eq. (^5]) we 
obtain 

i^{p) = \Y.^oki:s.k{p)-Xi{p) (32) 

kl 

Accordingly, we define the functional 

= ^ / E ^oki:^k (p) ■ X, (p) d^x (33) 

•' kl 
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Variations of ^'[p] imply variations of ^'[p, J] along paths defined by eq. ([5T|) . 
Now, the extremum property of ^[p] depends on the sign of 6^^^^ {p), that is, 

The same arguments presented for the problem (I^S)) can be used for eq. ([M]) . 
If /(?7) is positive definite it follows that S^^^'ipip) > so that ^[p] shows a 
minimum at the SNES. Thus the extremum properties of the functional vE'[p, J] 
along the two paths dealt with are determined uniquely by the sign of the form 
f{ri) which can be established only by thermodynamics. In this connection, we 
note that functional ([55)1 represent in a general form a Glansdorff-Prigogine local 
potential. Thus, if the extended MREP principle establishes the minimum of 
\E'[p] , automatically, it also establishes the maximum of ^'[p, J] along the fixed 
force path, thus allowing for an extended LDE principle. 

Let us consider the simplest case of a functional dependent on a single state 
variable. In the domain {X, J} we have -0 = J • X — (l/2)MoJ^ where Mq > 
because of the second thermodynamic law. The function describes a saddle- 
shaped surface where the saddle point {Xq, Jq} is a maximum along the path 
(5X — and a minimum along the path 6J — MqiSX. In the important case 
of the heat conduction we have X = V(l/r) and L = X{T)T'^ where A(T) 
is a temperature-dependent coefficient of thermal conductivity. The restricted 
functional is 

^[T]^\j KT,)T^[V{l/T)fd'x 
which is the local potential for heat conduction^. 



IV SOME EXAMPLES 



Bulk thermoelectricity in intrinsic semiconductors. We apply the above 
variational procedure to the case of bulk thermolectricity in semiconductors 
where cross effects between electron-hole and heat transport are to be taken 
into account. The entropy production is (compare with eq. [5] where — ^ is 
replaced by v) 

CTS = -Jn ■ V— -Jp • + Jq • V- + — (35) 

where and are the electrochemical potentials of electrons and holes, re- 
spectively, [tfn = p-n + e-V and ipp — fip — eV, V standing for electric potential, e 
for the absolute value of electron charge and /i for chemical potentials), v for the 
generation-recombination rate [v > for recombinations) and T = ipn + ipp for 
the carrier affinity (or, equivalently, for the differences of quasi- Fermi levels^). 
On these grounds, by taking into account of eq. p4)) (with fluxes as variables) 
we have 
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\MoqqJI + MonQ in ' Jq + 



T 



(Fx 



J„ • V — +jp 

By varying potential 4* with respect to temperature we obtain 



that is 



-Jn 



2^2 



2^2 



where Jq = Jq — ji so that (see eqs[5][6] ) 



(36) 



V • Js = crs 

As for the flux balance equations (variations with respect to ip), it is easy to 
verify that V • j„ = — u and V • jp = —v 

Now, let us consider the flux laws (variations with respect to fluxes). As for 
electron and hole fluxes, it is to be pointed out that we does not expect direct 
cross eff'ect. In this connection, later we will use (even if not explicitly shown) 



are 



, that is, MoqqMoup = Mq^qM^pq. The stationary fluxes 



Mopp T 



M, 



OnQ 



M, 



Onp . 



1 



1 



Monn 
MopQj 

Mopp 

MonQ 



MQnn 



Jn 



Mqpp 

MopQ 



(37) 
(38) 
(39) 



MoQQ ■ ^ MoQQ- MoQQ 

where the Onsager reciprocity relations will appear when the forces are explicitly 
shown. To write the flux equations in a more compact form it is convenient to 
define 



e 



= 1- 



M2 



MonnMoQQ 



MormAIl 



Ah 



M, 



OppMOQQ 



Opp 



(40) 



(41) 



Note that 1 - M^„q/Mo„„Moqq > and 1 - M^^q / MappM^QQ > due to 
positiveness of the quadratic form. Here a„ and ap stand for the electron and 
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hole conductivities, respectively. Thus, current densities of electrons and holes 
are respectively 

i.=j„ = _-(.v^-^ivr) (42, 
i = i, = -5(rvf-^ivT) ,43, 

which should be compared with the well known thermoelectric current equations^ 
^ = -^(v^n + f^VT) (44) 

^ = -5(v^. + §?Vr) (45) 

where Q„ and Qp stand for the excess of kinetic energy with respect to the mean 
energy of electrons and a holes, respectively. The following correspondence are 
found Qn = - {MonQ/MoQQ + fn) and Qp = - {Mqpq/Mqqq + ipp) .Based on 



these results, the heat flux equation can be reduced to the simple form 

jQ = -AVT + (Q„ + (^„)j„ + (Qp + (Pp)jp (46) 
where we used A = 1/MqqqT ^. Finally, the following relations are found 

Monn = (eVanT) + (1/A) [(g„ + iPn) /Tf (47) 

Mopp = {e'/apT) + (1/A) [{Qp + Vp) /Tf (48) 

- MonQ = {Qn + ^n) /A T ^ (49) 

- MopQ = {Qp + ^p) /A T 2 (50) 
Monp = {Qn + ) {Qp + ^p) /A T 2 (51) 



As for the local rate of electron-hole generation recombination process, we 
obtain v = M^J^T/T or, by putting Mq"^ = LowTq, v = LowT (at SNES). 



External forces 

The thermodynamics of SNES involve non isolated system because of the 
necessity of continuosly feeding the stationary processes. The external param- 
eters fixing the SNES are not included in the theory which only establishes the 
relations among physical quantities of the system dealt with. In several case 
this could generate difficulties if internal fluxes to be varied are bound with 
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the external ones. This problem can be figured out by means of proper ad- 
ditional constraining equations. However, we may also conveniently re-define 
the thermodynamic potentials to include explicitly the external fluxes. To this 
purpose, we use a different entropy balance equation in which negative terms 
are added to the entropy production. These additional terms account for the 
action of external forces continuously displacing the system from equilibrium. 
This is tantamount to say that we are considering the rate of the whole system 
entropy change, which is expected to vanish at SNES (the system entropy re- 
main constant). The nature of additional terms discriminate the internal fluxes 
(which appear in the kinetic term of potential -ijj) from the internal forces. In 
this connection, it is to be taken into account that derivatives with respect to 
forces Xa coniugated to external fluxes Za should satisfy 

- ^« (52) 



dXc, T 

Thus, we will consider as internal fluxes the ones which links to fixed ^a- Note 
that presence of external terms does not change the conclusions about the po- 
tential t/j which depends on the internal properties of the system. We have only 
a shift of the stationary point in the J, X map. Of course, this topic merges 
with the boundary problem discussed in section 5. 



The MelxwcU theorem of circuits 

To explain the above point, let us consider a network of resistances and 
voltage sources (thus a discrete system), for (example? ideal batteries.. All 
the chemical processes taking place inside a battery arc to be considered as 
external to the circuit (open system). The whole entropy of battery increases 
due to spontaneous processes, but it release a negative entropy —lada/T to the 
circuit due to the work of charge separation. The current crossing the voltage 
source depend on the circuit features. Thus, all the circuit currents are to be 
considered as forces in the potential ip. Instead, the corresponding affinities are 
to be considered as fluxes to be linked with fixed ^a- On this ground, potential 
ip should be written as 

^ = E^-E^-E^^o.r? (53) 

i a i 

SNES is to be searched for with additional current constraints, that is, 

^r = E^-E^-E \Mo.r'^ + E /MA.}, (54) 

i a i f3 

Alternatively, if we deal with current sources da (now the fixed external fluxes), 
SNES is to be searched for by means of 

^ = E^-E^-E^^o^^^ (55) 

i a i 
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and 



f^/-E^-E^-E ^^o.x? + ^2^1^ 5Mr{.}, r^„j) (56) 

i a i /3 

where gi3{T^iy,T^a}) a-re the affinity constraints. 

As example, let us consider the simple circuit formed by a battery ^ and 
two parallel resistances, Ri and i?2, connected to a resistance R3. Current 
constraints are I ^ Ii + I2 = Is- Thus, 

^^r = E^-^-E l^'^o^^' + Aa (/ - /2 - h) + Afc(/ - h) (57) 

i i 

From dflr/dT, = 0, dVtr/dh = and dVLr/dl it follows T, ^ h/ {T Mq,), 
Fi = = T Aa, T3 = T Xb and =ri + Tg = Ta + T3, respectively. As for 
potential ^ (for variations of along Sli = Moi STi paths) we have 

i 

where i?^ (T Moi)" • Equation ^ (or its generalization) by addition of 
current constraints call for the Maxwell theorem of circuits^^. 



Photo-excitations in semiconductors and recombination resistances. 

Now let us consider the electron-hole recombinations in semiconductors un- 
der steady injection of carrier in conduction and valence bands. For simplicity, 
we will deal with the case of uniform photo-excitation and temperature. Ac- 
cording to these settings we are assuming the infinite thermal conductivity of 
the lattice to which the heat produced by energy conversion (which accounts for 
the external entropy production) is released. Blackbody field is taken into ac- 
count in rate equations. Of course, the uniformity conditions make the foregoing 
considerations suitable for the system near equilibrium. Potential becomes 

= ^ - ^ - Imo^.v^ (59) 

where $ stands for the rate of (externally induced) carrier generation. The 
sign minus account for the negative contribution to the entropy induced in the 
system by the external force. Really, it is not relevant how the system is moved 
from equilibrium since we are dealing with a fixed rate $.^At the steady state 
we obtain the expected equation v = LoyyTo = $. As for tjj (for variations of ijj 
along 5v = Lqw^^ path) it holds 

~ iLo™r2 $r 

^=2^—-— (60) 
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Finally, it is easily recognizable that eq. ([55]) is similar to eq. ([55]) . Thus we can 
conclude that at SNES (near equilibrium) a Maxwell theorem can be called for 
the recombination rate distribution in a semiconductor analogously to the dis- 
tribution of currents in a electric ciurcuit; now resistances are to be replaced by 
recombination resistances related to the inverse of recombination probabilities^ 



V. FINAL REMARKS AND CONCLUSIONS. 

We still remark that the variational equation ([77)1 stands on the condition 
of having fixed state variables at the system boundary. We can remove this 
condition by considering a surface integral such that*. 



We need not modify anything in the conclusion of the previous section concern- 
ing the extremum properties of the functional 5*. Indeed, eq. (|6ip holds at any 
variation order so that the second variation of the modified functional coincides 
with (5(2)*. 

In the case grad v 7^ we must include the viscoelastic term in the entropy 
production. ^'^2 In this connection we point out that v cannot be considered as 
a functional variable. This is because the variational criterion (|27p is closely 
related to the entropy conservation law and does not account for momentum 
conservation. Thus, the functional (|33p. as it stands, cannot be used in the field 
of hydrodynamics. We need not go further in this matter which is beyond the 
scope of this paper. 

We have shown that in continuous systems satisfying LTE conditions, the 
(extended) MREP and (local) LDE principles are connected by means of a 
suitable functional. In the case of a single state variable, the density of the 
functional describes a saddle-shaped surface in the force-flux domain. The sad- 
dle point represents the stationary state. The MREP path is the one where 
the phenomenological coefficient maintains its reference state (SNES) value. In- 
stead, the LDE path is the one where the force is kept fixed. In particular, 
along the MREP path the functional behaves as a Glansdorff-Prigogine local 
potential. 
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